The Kohn mode for trapped Bose gases within the dielectric formalism 
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The presence of undamped harmonic center of mass oscillations of a weakly interacting Bose gas 
in a harmonic trap is demonstrated within the dielectric formalism for a previously introduced 
finite temperature approximation including exchange. The consistency of the approximation with 
the Kohn theorem is thereby demonstrated. The Kohn modes are found explicitly, generalizing an 
earlier zero-temperature result found in the literature. It is shown how the Kohn mode disappears 
from the single-particle spectrum, while remaining in the density oscillation spectrum, when the 
temperature increases from below to above the condensation temperature. 
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I. INTRODUCTION 

In 1961 W. Kohn (T) first noticed that the electron- 
electron interaction does not change the cyclotron reso- 
nance frequency in a bulk three dimensional gas. Brey et 
al. and Dobson |jj expanded the proof as a generalized 
Kohn theorem for parabolic quantum wells. Independent 
of the interactions the bare harmonic-oscillator frequency 
is preserved in the excitation spectra. These oscillations 
describe the center of mass movements along the main 
axes of the trap. They are used for calibration in the 
experiments with trapped Bose gases j|. 
Many different model approximations were derived since 
the discovery of Bose condensation in trapped dilute 
Bose gases. For them the satisfaction of the Kohn theo- 
rem provides an important consistency check. That the 
Kohn theorem is satisfied within the zero temperature 
Bogoliubov equations for a trapped dilute Bose gas was 
proven by Fetter and Rokhsar H. The Kohn theorem 
was demonstrated by Griffin and coworkers for various 
approximations leading to hydrodynamic equations Q. 
Stoof and coworkers checked the presence of the Kohn 
modes in an approximate variational approach to the ki- 
netic theory also in the collisionless regime @. A finite 
temperature theory based on a linear response formalism 
was given by Minguzzi and Tosi Q and the validity of the 
Kohn theorem in this approach was mentioned in their 
paper, but without proof. 

The purpose of the present paper is to examine the valid- 
ity of the Kohn theorem in a specific finite temperature 
approximation within the dielectric formalism. It was 
introduced and studied in a previous paper |J. It was 
shown there that its results agree with those results which 
can also be calculated in the Minguzzi and Tosi approach 
||, but that it gives additional results, in particular for 
the 1-particle Green's functions. 

The starting point is the closed set of eigenvalue equa- 
tions for the single particle Green's functions and the 



density-density correlation functions. We first identify 
the eigenmodes for the center of mass movement. Then 
we show that they are eigenfunctions of the density au- 
tocorrelation function at the bare trap frequency. 
We furthermore prove the validity of the Kohn theorem 
in the single particle spectra in the Bose condensed sys- 
tem and show how the Kohn modes disappear from the 
single-particle spectrum, but remain in the density au- 
tocorrelation spectrum if the number of particles in the 
condensate goes to zero. 



II. EIGENVALUE EQUATIONS FOR THE 
DENSITY AND SINGLE PARTICLE 
AUTOCORRELATION FUNCTIONS 

In this section we consider the eigenvalue equations 
for the single particle excitations and the density auto- 
correlation function above and below the critical temper- 
ature for Bose condensation T c . The eigenfrequencies ujk 
and the eigenmodes Sn to t,k(r) of the density correlation 
function \(r,T;r' ,t') = -~(T T [n(r,T)n(r l ,T')]), where 
h(r) = h(r) — (n(r)), are given by 



d 3 r' X - 1 (r,r',L0 k )6n tot>k (r') = 0. 



(1) 



It is useful to express \ m terms of its interaction- 
irreducible part x vm 

X{r,r',u) 

= X{r, r', <*>) + ! J d3r i X(r, n, w)x(ri, r', w) . (2) 

Here and in the following the tiWe-sign denotes 
interaction-irreducible quantities, which are given by di- 
agrams which cannot be split into two by cutting a single 
interaction line. 
Then Eq. dll) takes the form 
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Sn t ot,k(r) = t I d 3 r'x(r,r' ,uj k )5n to t,k(r') . 



(3) 



ft . 

Besides we consider the matrix of Green's functions 

G a ^(r, t; r\ r') = (T T [<p a (r, r)<p\{r' ', r')] ) (4) 

with field operators in Matsubara representation 
<pi(r, t) = <p(r,r) and <fii{r,T) = 0'(r,r). The cor- 
responding eigenfrequencies uj' k and eigenfunctions tp a ,k 
are determined by 



M-;(r,r>^ M M=o. 



(5) 



Of course we need to make approximations for x an( i 
G a p and describe now the model we previously intro- 
duced within the framework of the dielectric formalism 
P,fio| . The approximations will be further detailed in the 
subsections [I A and [I B and the appendix where they 
will be used to derive the eigenvalue equations for the 
different excitation spectra. At this point it is sufficient 
to fix the approximations for the thermal density m(r). 
The external trap potential U{r) is 



U{r) = \m J2 ^i x l 



(6) 



=1,2,3 



where are the main trap frequencies. 
The effective interaction potential v(r,r') has the form 
v(r,r') — gS(r — r') with the constant interaction 
strength g. It is taken to contain already the complete se- 
ries of ladder-diagrams describing repeated two-particle 
scattering, i.e. it is the 2-particle T-matrix. 
The stationary thermal density riT^r) is calculated using 
as auxiliary quantity the stationary Hartree-Fock Greek's 
function G HF 



Tioj - H Ht (r) G (r, r', u) = hS(r - r') . 



The Hartree-Fock Hamiltonian is 



# W = -t^A + U(r) — (i + 2gn tot (r) 
Zm 



(7) 



(8) 



with the chemical potential /i and the total density 
^tot(^) = nr(r)+n c (r), where n c is the condensate den- 
sity and is 



n T (r) 



]imJ2e iu)nr >G HF (r,r,uj n ) 



k B T 
Ti v 

n 



(9) 



A. The uncondensed Bose gas 

We first consider the uncondensed Bose gas. Here, the 
lowest order contribution to x is the bubble diagram x° 
containing two Hartree-Fock Green's functions 

X°(r,r» = 

~-k J2 GffF ( r ' r '' ^n)G HF (r\ r, iu> n - d) (10) 

n 

{n{£k) - n{£l) \ H k F {r')^ F (r') 



k.i 

i«-»k 



k.l 



1 



uj+ {e k -ei)/h 
1 



-lu + {e k - ei)/h 



(11) 



Here we first performed the Matsubara sum and then 
symmetrized the obtained expression with respect to 
k <-> I. 

X is represented, within the model, by summing up the 
multiple particle-hole scattering processes including ex- 
change with the result, in symbolic notation 



X(r,r',w) 



d i r 1 



l-|x»] (r.rO^^r'.w). (12) 



Inserting these diagrams x i n Eq. multiplying from 
the left with the integral operator [l — (g/~h)x°{u)] and 
combining similar terms, the density correlation spectra 
are given by 

*ntot,*(r) - 2| J d 3 r' X °{r,r',u; k )5n tot , k {r') . (13) 

Above T c the single particle Green's function in our 
model is just the Hartree-Fock approximation. Of course 
the corresponding eigenfrequencies oj' k = e k /h given by 

d 3 r'{G HF }-\r,r',u ) ' k )^ F (r')=^ F (r) (14) 

differ from those of the density correlation spectra u k . 

B. The Bose condensed gas 



u) n are the Matsubara frequencies 2mrkBT/h and n(e k ) 
is the Bose factor to the eigenenergy e k for the eigenfunc- 
tion (p k F (r) of the Hartree-Fock Hamiltonian. The latter 
will be taken as real, in the following, without restriction 
of generality because the Hamiltonian (g) is real. Within 
the model Q all internal lines of diagrams are calculated 
with G HF . 



In the presence of the condensate the two kind of exci- 
tations (|l]) and (||) share the same frequencies. This fact 
is a consequence of the underlying spontaneous symme- 
try breaking and can be conveniently formulated by using 
the dielectric formalism. This formalism expresses all the 
relevant quantities in terms of proper irreducible quanti- 
ties (whose graphs remain connected even if one cuts an 
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arbitrary, single interaction line or propagator). For the 
Bose condensed gas we have to include interaction pro- 
cesses with the condensate atoms, involving many fur- 
ther diagrams within our approximation. In order to 
avoid a somewhat lengthy detour we refer the interested 
reader to the appendix, where we briefly summarize the 
approximations and the derivation of the following equa- 
tions. These approximations have been presented in de- 
tail in reference (T^j . The complete set of eq uatio ns for 
the eigenfunctions within our model are (Eq. (|A12| )) 

(Ho(r) + gn c {rfj <p ljk {r) + g<f>l(r)p 2tk {r) = 

hu> k <Pi,k(r) - 2g<$> Q (r)5n TM (r) , (15) 
(&o(r) + gn c (r)^j ip 2 , k {r) + g<S>* 2 (r)<p 1:k (r) = 

-huj k ip2,k( r ) - 2g®o(r)5n T , k (r) , (16) 
where we defined 

£nc lfc (r) = $5(r)^i,fc(r) + $o(r)cp 2 , k (r) (17) 

and 

£n T , fe (r) = 6n totik (r) ~ Sn c ^ k (r) . (18) 

Due to the presence of the condensate we have quasiparti- 
cle excitations with two-component vector fields tpk,a{r) 
(a = 1,2) as eigenfunctions. 

The condensate wave function < &o( r ) is calculated 
from the finite-temperature generalization of the Gross- 
Pitaevskii-equation 



H (r)$o(r)=0, 
with the Hamiltonian Ho(r) 



(19) 



Ho(r) = - — A + U(r) - M + g\%{r)\ 2 + 2gn T {r) . 
Ira 

(20) 

We allow complex solutions of Eq. (p0|), i.e. our consid- 
erations apply also to the case where vortices are present. 
The total density is given by the sum tit + n c with (g) 
and n c (r) = |<I>o(r)| 2 . 

In the appendix we show (see Eq. (|A1C| )) that in the 
presence of a Bose condensate Eq. (|13[) must be rewrit- 
ten as 



SriT,k(r) 



,9 



d 3 r'x°(r,r',uj k )6n tot , k (r'), (21) 



which can be also expressed in terms of Sn Ctk 
5n T ,fc(r) = 2| Jd 3 ri Jd 3 r 2 



1 - 2f x°K) 



{r,r 1 )x°(r 1 ,r 2 ,u; k )6n c , k {r 2 ). (22) 



With Eq. @ we see that Eq. 



(jlf]) and <\l§ still form a 
set of homogeneous coupled equations for ip a , k . 



For n c = Eq. ( |l5| ) agrees with the Hartree-Fock equa- 
tion. In the low temperature region T — > we have 
X° — > and according to Eq. (pl|) Sut — > 0. Therefore, 
at T w Eqs. (|l|) and (|l|) approach the usual Hartree- 
Fock Bogoliubov Popov equations. 

Eqs. (|l|) and (H) agree with the linearization of the time 
dependent Gross-Pitaevskii equation around its station- 
ary solution $o(t) if the thermal density is assumed to 
be time-dependent and is also linearized around its sta- 
tionary value riT{r). Of course we then need additional 
equations like Eqs. (EEJ) and (Il7j) to fix <5nT,/c- 



III. THE KOHN THEOREM 

The Kohn theorem states that the center of mass os- 
cillations of the total density n to t along the main axes of 
the trap are unchanged by the interactions and appear 
at the bare trap frequencies £li as exact eigenvalues. 
For the proof we choose arbitrarily one of those directions 
ej. The center of mass movement is just a displacement 
of the total density in the trap. 

Therefore, the eigenfunctions to the Kohn frequency tti 
should be given by the infinitesimal displacement oper- 
ations Tjdi acting on the stationary densities (r) <C 1) 
r\ din{r) = n(r + r\ e^) — n(r) — 5n,i(r). 
The 

Kohn mode ^ntot,^'") should therefore be proportional 
to d l n tot {r) = din c (r) + 2 Y, k n(s k )<p k F (r)dip% F (r) if 
we insert the expression for n^. 

We prove that the Kohn theorem is satisfied in two steps 
where we set uo = fli in all the expressions. First we 
show that we get Sut^ = c^nx on the left hand side of 
Eq. ( pl|) if we insert 5ntot,i = ftritot on the right hand 
side. Using the result Stit^ — c^nx in Eqs. ( |i~5| ) and ( |l6| ) 
we calculate the corresponding Lp a ^ k from which we derive 
5n c .i — din c . We hereby reobtain 5n to t,i = diUtot proving 
that diUtot and (fa :k are eigenfunctions of the equations 
to the trap frequency Slj. In the following the effect of 
the partial differential operator on the eigenfunctions of 
H HF are evaluated by using the commutator relations 



0, , H HF (r ) = mtiixi + 2gd t n tot (r) , 



Xi,H HF (r 



-0, 



(23) 
(24) 



We obtain the matrix elements of the coordinate and the 
differential operators within the basis of the Hartree-Fock 
eigenfunctions 

(e k -e l )(l\x i \k) = —(l\d i \k), (25) 
m 

(e k -ei){l\di\k)=mn?(l\xi\k) 

+2g{l\d intot {r)\k) , (26) 

where we used the abbreviation 
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(27) 



In the first step we have to insert <5ntot,fc = c^ritot in the 
expression (|ll]) for the right hand side of Eq. (|2l]) with 

LU = Qi 

Sn T 4r) = (n(e k ) - n( Sl )) 



k,l 



(k\Sn totji \l) 



(k\5n to t,i\l) 

Mlj 4- (e k - ei) -Tifti + (e k - £/) 



(28) 



We further concentrate on the coefficients corresponding 
to^(r)^(r)(n(eO-nfo)) 



fJ 



(l\Sn tot< i\k) 



(l\Sn to t,i\k) 



= 2g 

(26) 



Mli + (e k - ei) -hQi + (e k - ei) 
(£k ~ £i)(l\Sn tot ,i\k) 



. (e k - e;) 2 - (^) 2 
(e k - £i) 2 (l\di\k) - (e k - £i)m^(l\xi\k) 



(e k - eiY - (hSlif 
(e k ~ ei) 2 (l\di\k) - (MU) 2 H\di\k) 



(e k - £l ) 2 - (nn t y 



= (l\di\k). 
We obtain the result 



(29) 
(30) 

(31) 
(32) 



Sn Tyi (r) 

= Y,m\k)^ F (r)^ F (r) (n(e k ) - n( £ ,)) (33) 



fe,i 



2j2n(s k )^ F (r)rf F (r)(l\d l \k), 



(34) 



fc,i 



where in the last step we used the antisymmetry of 
(l\di\k) under exchange of k and I. The expression (|34|) 
for 5nT t i is exactly the partial differential of nx with 
respect to Xi. Therefore, we get by comparing the left 
and the right hand side of Eq. ( pT| ) Sutj = din^. 
In the second step we have to insert Sn^^ — ftriT in 
Eqs. ([[5]) and ([Tq ). But then it is straight forward to 
see that Eqs. (|1E|) and (|T^) are solved for lo = Oj and 
5nT,i = diUT by setting 



Pi,i(r) = (di - XimQi/h)^ (r) , 
Va.i^) = (di + x l mVL i /h)<&l(r) . 



(35) 
(36) 



Inserting these expressions in Eq. (|I7| ) we get Sn c (r) — 
din c (r) and hence Sn to t(r) = Sn c (r) + Sti^fr) 
= dintot(v), completing the proof. Therefore the expres- 
sions dintot(r) and ip a .i(r) are eigenfunctions with the 
main trap frequencies f2i as eigenfrequencies. 
Fetter and Rokhsar g] proved that ([35]), (|3^) are the 
Kohn modes of the Bogoliubov equations at T = 0. Our 
preceding considerations generalize their proof to finite 
T. In the whole calculation we have chosen to fix the 



normalization of the density fluctuations Sntot,i = ft^tot- 
By performing a simple partial integration we get the 
normalization of the Kohn eigenfunctions fi^(r), f2,i(r) 
according to / d?r(\cpi ti (r)\ 2 - \^2,i(r)\ 2 ) = 2mfliN c /K, 
where N c is the total number of atoms in the condensate. 
With N c — ► the Kohn modes are seen to disappear from 
the single-particle spectrum. 

For N c — (i.e. $ = = n c ) Eq. (JTsJ) agrees with 
the usual Hartree-Fock equation where the Kohn mode 
is no longer present in the single particle spectra. But 
the Kohn modes survive in the density correlation spec- 
tra given by Eq. ( ]2l| ) for rj to t = ^t- In this case the 
proof of the Kohn theorem is already finished after the 
first step of our proof. 



IV. DISCUSSION AND SUMMARY 

The fundamental difference between a spatially homo- 
geneous and a trapped Bose gas is the absence of the 
conservation las of momentum in the latter, which plays 
an important role in the theory of homogeneous systems. 
Certainly any approximation made in homogeneous sys- 
tems should satisfy this fundamental law. In the special 
case of harmonic trapping potentials the momentum con- 
servation law is replaced by Kohn's theorem, which states 
that there are three exact special modes where the cen- 
ter of mass oscillates harmonically with the three main 
trap frequencies. While this is easily proven for the ex- 
act Hamiltonian, the question whether any given approx- 
imation still satisfies this theorem is, in general, a rather 
nontrivial one. In the case of Bose condensed gases the 
question acquires an additional aspect, because of the co- 
incidence of single particle and density oscillation spectra 
in Bose condensates. The Kohn mode, which is a density 
oscillation mode, must then have a single particle coun- 
terpart, which somehow disappears from the spectrum if 
the temperature is raised, so that the condensate disap- 
pears and the single particle spectrum and the density 
oscillation spectrum become decoupled. 
In the present paper we have examined these questions 
for a specific approximation which includes direct inter- 
action and exchange and is formulated within the dielec- 
tric formalism, which guarantees form the onset the co- 
incidence of single particle and density oscillation spec- 
tra in the Bose condensed regime. By deriving a closed 
set of equations (|llj)-([T^) for the single particle modes 
and the density modes generalizing the usual Bogoliubov- 
deGennes equations in Popov approximation, and solv- 
ing them for the special coupled modes, corresponding 
for the densities just to a translation, we could not only 
verify the Kohn theorem for the closed set of equations, 
but obtain explicit expressions for the single particle 
modes ([35]), (pf). These expressions for the latter are 
beautifully simple. They generalize the corresponding 
zero-temperature result obtained by Fetter and Rokhsar 
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by replacing the zero-temperature solution of the Gross- 
Pitaevskii equation by its finite counterpart of our spe- 
cific model within the dielectric formalism. In this form 
our results (|35|), ( |36|) show explicitly how the single par- 
ticle component of the Kohn mode vanishes while keep- 
ing the density oscillations component unchanged, if the 
number of particles in the condensate is sent to zero. 
A simpler version of the model treated here leaves out 
all exchange processes. This is then the simplest model 
of an interacting Bose gas which one can set up within 
the dielectric formalism. It was analyzed in detail for a 
homogeneous Bose gas in 11 1, but one can also analyze 
it for the trapped gas along the lines of || extending 
work in Using the same procedure as employed in 
the present paper one can readily show that the Kohn 
theorem is respected also by this simpler model. 



In the model introduced in ||,[T(| the vertex functions 
Aq 1 ^ contain the lowest order contributions given by 
$* (r) and a series of higher order diagrams correspond- 
ing to the diagrams summed up in Eq. ( |l2| ) 

AW(r,r', W ) = $* )0 (r)5(r-r / ) + | f d 3 ri A^(r, ri ,Lj) 



xx (ri,r',u)) 



(A4) 



If we allow $o to be complex we have to use $o for the ex- 
citation processes due to the corresponding annihilation 
of a condensate atom and $q m the opposite case. In 
order to simplify the notations we write <& Q .o ( a = 1, 2) 
with $i i0 = $0 and <& 2 ,o = 

~ i'f) 

G a f3 is in general expressed by the regular part of 
the self-energies in the form 



APPENDIX: 

In the first part of the appendix we introduce the nec- 
essary quantities given by our model PJlol including ex- 
change to describe the density autocorrelation function 
for T < T c . Then we show how Eqs. @ , (|l6|) , which 
we have used here as a convenient starting point, are ob- 
tained. We first note that the contributions to x above T c 
given in Eq.(p"2|) are not only interaction line irreducible 
but also propagator line irreducible. We denote these 
so-called reguiar contributions by x^ ■ Below T c we get 
contributions to x i n addition to the regular ones, namely 

X (r,r',uj) = x (r) (r,r',w) + J d 3 r ± J d 3 r 2 A ( a r) (r,r u uj) 

xG^(r 1 ,r 2 , W )Aj ) *(r 2 ,r', W ). (Al) 

The interaction line irreducible Green's functions G a fj 
are related to the exact Green's functions G a p by the 
Dyson-type equation 

G a ,0{r,r',u;) = G a .p{r,r',uj) 

+ f/d 3 ri J d 3 r 2 J d 3 r 3 G a , 1 (r,r 1 ,iu) 

xA*(ri, r 2 , cj)A 5 (r 2 , r 3 , u;)G St0 (r 3 , r', u) , (A2) 

where the vertex functions A a are related to their regular 
parts by 

A a (r,r',u) =A a (r,r',uj) 

+ §/ d 3 ra W (r,r 1)W )A a (r 1 ,r', U ). (A3) 

The regular vertex functions A& express the coupling of 
the single-particle Green's functions G a /3 to the density 
autocorrelation function x- The coupling comes from ex- 
citation processes out of and relaxation processes into the 
condensate. 



/ huj 



(r, r') — S(r — r') 



2rn +V-U{V) 









-Hlo 



2m +»-U(r)J 
(A5) 



The approximation of the model for £( r ' is defined by 
E<g(r, r', u) = 6(r - r')(gn c (r') + 2gn T (r'))S a0 , 



-^$ Q ,o(r)x (r) (r,r',c)^ (r'). 



(A6) 



Like in Eqs. (|12[),( A4) we add to the lowest order regu- 
lar self-energy diagrams gn c (r') + 2gnr(r') appearing in 
the diagonal of the sum over multiple particle-hole 
scattering processes including exchange with additional 
factors <f>5 1 due to excitation processes out of and the 
absorption processes into the condensate. 
Below T c the eigenf unctions fi(r) and ip2(f) can be ex- 
pressed as a linear functional of the corresponding eigen- 
function Sn t ot.k( r ) of the density autocorrelation func- 
tion spectra (see |12|) 



ip a ,k( r )= J <^ r i J d 3 r 2 ^2G a p{r,r 



xA {r} *(r 1 ,r 2 , w fc )5n t ot,fe(r 2 ) . 



(A7) 



Without any approximations for the regular quantities 
S^g, x^ an d Aa ' the functions (p a {r) would fulfill Eq. 
® if <5n totlfe (r) fulfills Eq. (§. 

We insert Eq. (Al) in Eq. (|J) and take into account 
~1) to get 



Sn to t,k(r) 



^X {r) (r,r',uj k )Sn tot:k (r') 



d 3 r' 



(A8) 
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Defining 5n c , k (r) = %(r)<p ltk (r)+$ (r)(p 2 , k (r) and sep- 

a (r) 

arating out the the lowest order diagram of A a we 
derive an equation for SriT,k = ^tot,fe — Sn Cj k 



fe T , fe (r) = | y d 3 r'x (r) (r,r',u k ) 

x ((Jn c> fc(r') + <5n t ot : fc(r')) 



(A9) 



Eliminating Sn c , k by using #n C) fc = <5n to t,fc — <5ft"r,fc this 
equation is solved by 

*nr,*(r) = 2^ f d 3 r'x°(r, r', ^)5n tot , fc (r') , (AlO) 



as can be verified by using the result ( pL 2| ) for x ■ 
Acting with the operator j-G(uj k ) (r,r') on tpp tk (r') 

a/3 



and using (Al)-(A7) we obtain 
{±Tiw k - H (r))(pi fe (r) 



efV$i (r)x (r) (r,r',a; fc )<yTi c , fc (r') 



+^i i0 (r) y dV (5(r-r') + |x (r) (^,r-',c)) 
x<5n totifc (r') , (All) 



with the help of Eq. 
results 



where Hq is given by ( 20h. Eqs. (All) can be simplified 



A9|) and we can directly give the 



{±huj k - H Q (r))ipi_ k (r) = 
g$i (r) (5n to t,fc(r) + 5n T ,fc(r)) . 



(A12) 



Eqs. (A12) agree with Eqs. (|T^) and ( p^ ) using the 
definition 8n c , k = $oVi,fc + 
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